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Abstract. 

We present the results from a suite of N-body simulations of the tidal stripping of two- 
component dwarf galaxies comprising some stars and dark matter. We show that recent kine- 
matic data from the local group dwarf spheroidal (dSph) galaxies suggests that dSph galaxies 
must be sufficiently massive (10 9 — 10 10 Mq) that tidal stripping is of little importance for the 
stars. We discuss the implications of these massive dSph galaxies for cosmology and galaxy 
formation. 
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In the current cosmological p aradigm all structure forms from the successive mergers 
of smaller substructures (see e.g. lWhite fc Reesll 19781) . Cosmological N-body simulations 
demonstrate that much of this substructure w ill survive up to the present epoch in 
i-^ ' the form of dwarf galaxies ijNavarro et aljfl996() . While it is well known that cold dark 
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1. Introduction 



matter (CDM) theories over-produce this substructure, this 'sub-structure problem' can 



be solved by either suggesting that only the most massive s ubstructure halos form s tars 
£j ! JStoehr et all 120021): or by suppressing small scale power l|Avila- Reese et, aT1l200ll and 

IZentner fc BulloclJl20ol . 

Either way, these solutions suggest that dSph galaxies in the local group must be 
massive - some ~ 10 9 — 10 10 Mo in total mass. At this mass, for all but the most extreme 
, plung ing orbits, tidal stripping would not affect the central stars and gas ijHavashi et alJ 

l2ooa> . 

In contrast to this picture, there is compelling evidence to suggest that tidal strip- 
ping could be extremely important for t he formation and evolut ion of the local group 
dSph galaxies. IWalcher et all l|2003f) and IWilkinson et all l|2004l) point out that all of 
the local group dSph galaxies observed with deep enough photometry have a break in 
the light profile which i s similar to that obser ved in tidal stripping simulations (see e.g. 
Aguilar fc White! Il986t Ijohnston et alJl2002l and iGomez-Flechoso fc Martmez-Dekadol 



2003|h Furthermore. iMaver et alJ <|2001bj) and lMaver et all l|2001a() have suggested that 



the very low internal angular momenta and spheroidal morphologies observed in local 
group dSph galaxies could be explained by tidal effects. However, to experience strong 
tidal effects, the local group dSph galaxies must reside in low mass dark matter halos 
( ^ 10 9 M(n) which appears to b e at odds with naive co smological predictions. 

iKazantzidis et all l)2004() and iKravtsov et all l)2004|) have rece ntly argued tha t dSph 
galaxies might not inhabit the most massive substructure halos. IKravtsov et alJ l)2004j) 
suggest that reionisation sets an epoch at which the smallest collapsed structures cannot 
form stars, rather than having a mass scale below which star formation ceases. This can go 
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so me way towards s olving the mystery, but not all of the way: even in the scheme proposed 
bv lKravtsov et all 1 2004^ dSphs must inhabit the more massive halos on average. 

In order to resolve these issues we r eally require more obse rvational constraints from 
the local group dSph galaxies. Recentlv lWilkinson et alJ ((2.004) have found a sharp fall-off 
at ~ 1 kpc in the projected velocity dispersion of both the Draco and Ursa Minor (UMi) 
dSph galaxies. At first glance, these data would seem to rule out the tidal stripping 
model. It is w ell known that tid al stripping produces heating of stars within the tidal 
radius (see e.g. lOh et alJIlQQ^ and lcTomez-Flechoso fc Martme~ Dfilgadffennaji . Thus the 
existence of a cold population at ~ 1 kpc suggests that the tidal radius must lie beyond 
this point and that dSph are then 'high mass' galaxies. 

In this paper we examine this argument more carefully. We present the results from 
a suite of N-body numerical simulations of the tidal stripping of two component dwarf 
galaxies comprising some stars and dark matter, in orbit around a Milky Way-like galaxy. 



2. Results 

The initial conditions for the dwarf galaxies were set up as in lKazantzidis et alJ l|2004f) . 
but for two-component spherical galaxies comprising a dark matter halo and some stars. 
The satellite galaxies were well resolved with 10 5 star and 10 6 dark matter particles. 

The N-body integration for the orbit of the dwarf galaxy about the Milky Way was 
then performed using a modified version of the GADGET N-body code ^Springe! et alJ 
l200lh . The code allowed for a fixed potential to be introduced to model the Milky Way. 
We used a Miyamoto-Nagai potential for the Milky Way disc and bulge and a logarithmic 
potential for the Milky Way dark matter halo. The simulations, in itial conditions and 
model parameters are described in more detail in Ifiead et al. (200^). 

Figure ^ shows the results for three very different dwarf galaxies on orbits around the 
Milky Way after ~ 10 Gyrs. In each case the orbiting galaxy is currently at apocentre. 

The first model (model A) is one in which the satellite g alaxy contains no dark mat- 
ter. It was a Plummer sphere l)Binnev fc Tremainel ll987T with total mass and scale 
length of 10 7 M Q and 0.23 kpc. The galaxy was placed on a plunging orbit with apoc- 
entre and pericentre of 80 and 20 kpc respectively. Notice, as published by previous 
authors, that the isodensity contours (left panel) show the tidal 'S-shape' caused by 
the leading and trailing tidal debris. Several authors have argued that there is evidence 
for such shapes in the isodensity contour s of the local group d Sph galaxies (see e.g. 
iGomez-Flechoso fc Martinez- Delgadol2003l and iMnnoz etahl EoOM : and that this, there- 
fore, implies that the local group dSph galaxies must be strongly tidally distorted. Con- 
sider, however, the middle panel for the line of sight velocities. There is a strong velocity 
gradient across the th e galaxy which is not observed in Draco or Ursa Minor (UMi) 
^Wilkinson et al.ll2004f) . Furthermore, the projected velocity dispersion is rising beyond 
the tidal radius, where as in both Draco and Ursa Minor it is observed to drop sharply 
l)Wilkinson et al.|l2004h . Projections which give distorted isodensity contour s, small ve- 
locity gradients and near-flat projected velocity dispersions are very rare ijRead et alJ 
l2005(l . In no projection was the projected velocity dispersion of the stars found to fall. 

Even if we are viewing the local group dSph galaxies from a projection which minimises 
any kinematic trace of tidal stripping, notice that the central velocity dispersion in model 
A is much smaller than th at observed in the D raco dSph galaxy (see data points with 
large crosses). As shown bv lKlevna et alJ l|200lf) . the only way to reproduce such a large 
central velocity dispersion with the s tars we see in Drac o is to introduce a significant 
dark matter component (and see also lKlessen et aJj |2003j). 

Model B contained a significant Plummer dark matter component with mass 14 x 
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Figure 1. Results from three models: A (top); B (middle) and C (bottom). Left panels show 
the projected surface density of the stars, the middle panels show the line of sight velocities of 
the stars and the right panels show the projected velocity dispersion of the stars as a function 
of radius. The contour bars beneath the top left and middle panels show the number of star 
particles and the projected velocity in km/s respectively. Over-plotted on the righ t panels (large 
crosses) are the data points from the Draco dSph galaxy (Wilkinson et al. 20041). 

10 7 M Q and scale length 0.5 kpc; its orbit was identical to model A. Even after lOGyrs 
the galaxy has regular isodensity contours, no velocity gradient and a falling projected 
velocity dispersion. This is because the initial mass and concentration of dark matter was 
chosen to be sufficiently high that tidal stripping produced almost no effect on the stars 
over 10 Gyrs. Notice that for this model, the central velocity dispersion is much closer to 
that of Draco. 

Finally, model C was chosen to have a dark matter halo consistent with those observed 
in cosmological simulations. We used a Hernquist profile with mass 1O 1O M0 and scale 
length 1.956 kpc. In order for it to have a final velocity dispersion low enough given this 
starting mass we had to choose an extreme orbit with a pericentre of just 6kpcj\ At 
this radius the galaxy would likely affect the Milky Way disc (and dynamical friction 



f The central velocity dispersion is much lower than in the initial conditions for this model 
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would become important) and so already the orbit is quite unrealistic. However, notice 
that even for such an extreme orbit, the central velocity dispersion is still too large after 
10 Gyrs to be consistent with Draco. 

3. Conclusions 

In conclusion, tidal stripping cannot be very strong for many, if not all, of the local 
group dSphs. Strong tidal stripping, which would produce distorted isodcnsity contours, 
also leads to velocity gradients and flat or rising projected velocity dispersions - neither of 
which are observed in the local group dSphs for which we have good kinematic data (but 
see also lMunoz et al.12 005'1 . This suggests that dSph galaxies must be sufficiently massive 
such that tidal stripping is of little importance for the stars. Either they are on orbits with 
large pericentres, in which case they can have masses as low as ~ 10 8 M Q (jKlevna et alJ 
2001); or they are on more extreme orbits in which case they must be ~ 10 9 — 1O 1O M0 
depending on the extremity of the orbit. Our current cosmological paradigm would favour 
the latter hypothesis, but this leaves us with a puzzle: if the dSph are really as massive 
as ~ 1O 1O M0 and have dark matter densities which are cosmologically consistent then 
they would have central velocity dispersions which are too large to be consistent with 
Draco or UMi - even after significant tidal stripping and shocking. 
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due to disc shocking and tidal stripping (see Read et al. 2005 for more details). It does tend to 
zero at the centre as it should, but this is disguised somewhat by the choice of bin size. 



